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Abstract 

We consider strongly coupled gauge theory plasma with conserved global charges that 
allow for a dual gravitational description. We study the shear viscosity of the gauge 
theory plasma in the presence of chemical potentials for these charges. Using gauge 
theory/string theory correspondence we prove that at large 't Hooft coupling the ratio 
of the shear viscosity to the entropy density is universal. 
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1 Introduction 



The gauge theory/string theory correspondence of Maldacena [1,2] provides a valuable 
insight into the nonperturbative dynamics of strongly coupled gauge theory plasma. 
Originally, the string theory correspondence was formulated for static properties of 
strongly coupled gauge theories. It was pointed out in [3] that computing equilibrium 
two-point correlation functions of stress-energy tensor one can extract the physics of 
the near-equilibrium (hydrodynamic) description of strongly coupled hot gauge theory 
plasma. The computations of the finite temperature transport properties of the Af = 4 
SU(N) supersymmetric Yang-Mills (SYM) theory plasma at large t' Hooft coupling 1 
[3,6] were extended to various non-conformal gauge theory plasmas in [7-15]. One of 
the most surprising results of these computations was the discovery of the universality 
of the shear viscosity of strongly coupled gauge theory plasma. Specifically, for all 
gauge theories which allow for a dual supergravity description, and without chemical 
potentials for conserved global charges (if such charges are present), it was shown that 
the ratio of the shear viscosity rj to the entropy density s is a universal constant [16-18] 



in units H — ks — 1. 

Af = 4 SYM has a maximum t/(l) 3 C SO(Q)r abelian subgroup of the R-symmetry 

for which one can introduce (at most three) different chemical potentials. At strong 

coupling and finite temperature, the dual supergravity description of this gauge theory 

plasma is given by a system of near extremal D3 branes with generically three different 

angular momenta along the five-sphere S 5 [19] 2 . This supergravity solution allows for 

a consistent Kaluza-Klein reduction to five dimensions, where it is described within 

Af = 2 supergravity coupled to two abelian vectors commonly referred to as the STU 

model [21]. Within this effective five dimensional description, finite temperature Af = 4 

SYM plasma is dual to non-extremal black holes carrying generically three different 

£7(1) charges corresponding to three different S 5 angular momenta of the near-extremal 

D3 branes in the type IIB supergravity description. Even though Af = 4 SYM plasma 

with chemical potentials violates the assumptions of the universality theorem [16-18], 

explicit computation of shear viscosity leads to (1.1) [22-24]. 

1 Finite 't Hooft coupling corrections to Af = 4 hydrodynamics were discussed in [4,5]. 
2 The misprints in the five-form expression for this supergravity solution were fixed in [20]. 
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Expectation that (1.1) might in fact be more universal than anticipated in [16- 
18] was strengthened 3 by explicit construction of new models of finite temperature 
gauge/string theory correspondence with an _R-charge chemical potential [26]. Specif- 
ically, it was found that shear viscosity of all strongly coupled Y p ' q quiver gauge the- 
ory [27-30] plasmas with a U(l) R-charge chemical potential satisfies the universal 
relation (1.1). It was further conjectured in [26] that (1.1) holds true even in the 
presence of chemical potentials. 

In this paper we prove that (1.1) is indeed satisfied for generic strongly coupled 
gauge theory plasma with global charge chemical potentials. We begin with spelling 
out explicit assumptions (and their implications) under which we obtain (1.1). 

■ First, we consider strongly coupled gauge theory plasma that allow for a dual string 
theory description. This allows us to use the gauge theory/string theory Minkowski- 
space prescription [31] for computing two-point correlators of the stress-energy tensor. 
We will work in the regime of large 't Hooft coupling, where the supergravity approx- 
imation to string theory is valid. 

■ The gauge theory/ string theory correspondence relates a (_D+l)-dimensional strongly 
coupled gauge theory on M. ' 1 space-time to a particular background of ten-dimensional 
type IIB supergravity. In all known examples of gauge theory/string theory correspon- 
dence it is possible to do a Kaluza-Klein reduction along a compact (8 — D) -dimensional 
manifold and relate (D + l)-dimensional strongly coupled gauge theory to an effective 
(D + 2)-dimensional gauged supergravity, obtained from a consistent truncation of the 
full ten-dimensional supergravity. However, we do not believe that this is always pos- 
sible to do. In fact, the strongest form of the universality (1.1) proven in [18] does 
not rely on this assumption 4 . In this paper we assume that our gauge theory plasma 
is described by a type IIB gravitational background that allow for a consistent trun- 
cation to a (D + 2)-dimensional gauged supergravity. However, we will not assume 
that the geometry is asymptotically flat, or that the asymptotically flat region can be 
"re-attached". Thus, our universality class of (1.1) is necessarily weaker than that of 

3 Shear viscosity of the M2-brane plasma also appears to satisfy the universal relation (1.1) [25]. 
4 Both proofs [16, 17] implicitly assume the existence of a consistent truncation. The proof [17] 

further implicitly assumes that the background is asymptotically flat, as it relies on the universality of 
the cross section for the graviton scattering from the black hole horizon [32,33]. The latter universality 
was derived only for asymptotically flat space-times, and in general it is not known how to "re-attach" 
the asymptotically flat region in gravitational dual to non-conformal gauge theories (other than those 
described by flat Dp branes). 
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the zero chemical potential case considered in [18]. Since we would like to study gauge 
theory plasma with finite chemical potentials, these gauge theories must have a set of 
abelian conserved charges — otherwise we simply would not be able to introduce a 
chemical potential. In the dual supergravity picture, for each conserved U(l) charge, 
there must be a corresponding U(l) isometry. Introducing a chemical potential for a 
particular U(l) charge on a gauge theory side results in gauging corresponding U(l) 
isometry on the supergravity side. In the effective (D + 2) -dimensional gravitational 
description, besides Einstein-Hilbert term and scalar fields, one gets a set of Maxwell 
fields, one for each gauged isometry (or a chemical potential from the gauge theory 
perspective). To summarize, we consider strongly coupled (D + l)-dimensional gauge 
theory plasmas which have a following dual effective (D + 2)-dimensional gravitational 
description 

Sd + 2 = / dt D+2 V=gU - IC^d^dy? - r ab (<f>)F$F^ _ v(0)l 

16nG D+2 Jm d+ 2 I J 

(1.2) 

where V, 1C a p and r ab are arbitrary functions of scalar fields <p a , and the index (a) in 
Maxwell fields F$ = d^Aft - d u A ft runs over the set of nonzero chemical potentials. 
Corresponding to finite temperature (D + l)-dimensional gauge theory with chemical 
potentials, the effective action (1.2) must admit a black .D- dimensional brane solution, 
electrically charged under vector fields A a . The SO(D) invariance of the solution 
implies that it is of the form 

ds 2 D+2 = - c?(r) (dtf + cl(r) f>a;*) 2 + c^(r) (drf , 

i=i (1-3) 

Aft=5l^\r), Q = 0», 
for some scalar potentials 

■ Finally, we assume that the black brane horizon of (1.3) is nonsingular. This implies 
that there is a choice of the radial coordinate such that as r — > r ho 



wrizon 



ci -> «i(r - r horizon ) 1/2 , c 2 -> a 2 , c 3 -> a 3 (r - r^ ori20 „) 1/2 , (1.4) 
where constants satisfy 

«i «2 «3 7^ . (1.5) 
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Notice that given (1.4) the black brane temperature T and the entropy density s are 
correspondingly 



47ra 3 4G D+2 



2 The proof 



In this section, we examine the hydrodynamics of the gauge theory plasma at finite 
chemical potential dual to the generic black hole solution (1.3). In particular, using 
prescription [31], we compute the retarded Green's function of the boundary stress- 
energy tensor T^ u {t,x % ) (// = {t,x 1 }) at zero spatial momentum, and in the low-energy 
limit uj — > 0: 

Gf 2)12 (w,0) = -i J dtd D x e^9(t)([T 12 (t,x l ),T 12 (0,0)]) . (2.1) 

Computation of this Green's function allows for a determination of the shear viscosity 
T] through the Kubo relation 

V = lim ^ [Gft >12 (u;, 0) - G* >12 (a;, 0)] , (2.2) 

where the advanced Green's function is given by G a {uj,0) = (G R {u,0)y. We find 

<?£>;, 0) = -^(l + o(£)) , (2.3) 

where s is the entropy density. Inserting this expression into (2.2) yields the universal 
ratio (1.1). 

We begin the computation of (2.1) by recalling that the coupling between the 
boundary value of the graviton and the stress-energy tensor of a gauge theory is given 
by 5g\Tl/2. According to the gauge/gravity prescription, in order to compute the 
retarded thermal two-point function (2.1), we should add a small bulk perturbation 
5gi2(t,r) to the metric (1.3) 

ds 2 D+2 -»• ds D+2 + Sgi 2 (t, r) dx 1 dx 2 , (2.4) 

and compute the on-shell action as a functional of its boundary value 5g\ 2 it) ■ Symmetry 
arguments [34] guarantee that for a perturbation of this type in the background (1.3) 
all the other components of a generic perturbation Sg^, along with the gauge potentials 
perturbations 5 Aft and scalar perturbations 5<p a , can be consistently set to zero. 
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Instead of working directly with 5gi 2 , we find it convenient to introduce the field 
ip = ip{t,r) according to 

^ = 5 9i2 = 2 S 9i2 ■ (2.5) 

The retarded correlation function Gi 212 (cu,0) can be extracted from the (quadratic) 
boundary effective action Sb OU ndary for the metric fluctuations ip h given by 



d D+1 k dMn+2 
S boundary^] = / 4> b (-Uj) F(u , r) 4> b (lj) 



horizon 



where 



In particular, the Green's function is given simply by 



(2.6) 



(2.7) 



8M 



D+2 



G« tia (u>,0)= lim 2F(u>,r), (2.8) 

9M r D+2 -^dM D+2 

where T is the kernel in (2.6). The boundary metric functional is defined as 

S bounda r y m = HMr hm (s[ ul M + + 5— [fl) , (2.9) 

where Sl ulk is the bulk Minkowski-space effective supergravity action (1.2) on a cut- 
off space M. r £, +2 (where M.£> +2 in (1.3) is regularized by the compact manifold M. r D+2 
with a boundary dM. r D+2 ). Also, Sqh is the standard Gibbons- Hawking term over the 
regularized boundary dM. r D+2 . The regularized bulk action S[ u[k is evaluated on-shell 
for the bulk metric fluctuations ip(t,r) subject to the following boundary conditions: 

lir * V^,r) = /(t), 
9M r D+2 -*dM D +2 ^2 io) 

(b) : ip(t,r) is an incoming wave at the horizon. 

The purpose of the boundary counterterm S counter is to remove divergent (as dM r D+2 — > 
<9.Md+2) and ^-independent contributions from the kernel J 7 of (2.6). 

The effective bulk action for ip(t,r) is derived in Appendix A. It takes the form 



''bum 



_1 f d D+2 ti Cd+2 = f d D+2 ^ 

G D+2 J Mn ,„ ~ lQnG D+2 ju„,, 



l6nG D+2 J Md+2 Jm d+2 



Cl / V c 3 C 3 



The second line in (2.11) is the effective action for a minimally coupled scalar in the 
geometry (1.3), while the third line is a total derivative. Thus the bulk equation of 
motion for tp is that of a minimally coupled scalar in (1.3). Decomposing ip as 

^(t,r) = e"^W, (2.12) 
we find that the equation of motion reduces to 



r 2 



l n £i^L 

c 3 



-^ = 0, (2.13) 



C 3 



where primes denote derivatives with respect to r. 

Let's understand the characteristic indices of ip^ as r — > rhorizon- Consider the 
following ansatz 

^(r) ~ (A(r-r ho „, on )) A , (2.14) 

where A is a characteristic, finite, ^-independent energy scale associated with the 
background geometry (1.3). For background geometries dual to conformal gauge theory 
plasmas the role of A is being played by the temperature. More generally, A could be a 
scale at which conformal invariance is broken (either explicitly or spontaneously). It is 
crucial for A to be finite — otherwise the hydrodynamic approximation to a strongly 
coupled gauge theory plasma would not applicable 5 . Substituting (2.14) into (2.13) and 
using the near horizon asymptotics (1.4), we find to leading order in (A(r — rhorizon)) <C 
1, 

(r - r hortzon )^ K + ^^ + 4i=0, (2-15) 



2 

CVq Q^q 



or 



A = ±i = (2.16) 
where we used (1.6). Given (2.12), for the incoming wave at the horizon we must have 
^incoming „ (A(r _ r honzon) y^T) , (A(r - r honzon )) « 1 , (2.17) 

which to leading order in ^ takes the form 

^incoming „ j _ f " j n((A(r _ ( 2 .1 8 ) 



5 This applicability of the hydrodynamic regime is an implicit assumption of all the studies of 
strongly coupled gauge theory plasma within gauge theory/string theory correspondence. 
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valid when both inequalities are satisfied 



ln((A(r - r horizon )) 



> 1 



& 



to 
T 



ln((A(r - r horizon )) 



< 1. 



(2.19) 



Clearly, for finite A and sufficiently small uo there is an overlap region in (2.19). 

In what follows we will need a solution to (2.13), subject to (2.10), to linear order 
in ^, i.e., in the low-frequency approximation. In this approximation one can neglect 
the first term in (2.13) and write down the most general solution as 

Jr Ci(p)c 2 (p)^ 

where Ai{uj) are constants, depending at most linearly on ^. The first boundary 
condition in (2.10) implies that 

Ai(u) = l. (2.21) 
Substituting (1.4) into (2.20) we obtain the leading near horizon behavior of ^ 
a 3 A 2 (ou) 



1 - 



— ln(A(r - r horizon )) , 



ln(A(r -r horizon )) 



» 1. 



(2.22) 



Comparing (2.22) with (2.18) in the hydrodynamic approximation, i.e., for sufficiently 
small ^, we conclude that 



A 2 (u) 



i a® oj . 



To summarize, we have to leading order in ^ 

poo 

■0w( r ) — 1 + i a 2 ^ 

Notice that as r — > oo 



c 3 (p)rfp 
ci(p)c 2 (p) fl • 

c 3 (r) 



(2.23) 



(2.24) 



(2.25) 



c 1 {r-)c 2 (j-) D ' 

Once the bulk fluctuations are on-shell (i.e., satisfy equations of motion) the bulk 
gravitational Lagrangian becomes a total derivative. From (2.11) we find (without 
dropping any terms) 

Cd+2 = d t J l + d r J r , (2.26) 



where 



J = 



3c 2 c 3 _ 



2ci 



J' 



3c2 Cx 
2c, 



Cic: 



(2.27) 



2 _/.2 



C 3 



Additionally, the Gibbons-Hawking term provides an extra contribution so that 

jr _ jr _ ^1 ^ _ (^47 (2>2g) 
C 3 c 3 

We are now ready to extract the kernel JF of (2.6). The regularized boundary 
effective action for tp is 

SboundaryW = T7~~ri / dtd°X ( ^d r ^] + C.t. , (2.29) 

lb7lG D+2 J 9 M- D+2 \ 2 °3 J 

where, as prescribed in [31], we need only to keep the boundary contribution. In (2.29) 
c.t. stands for (finite) contact terms that will not be important for computations. 
Substituting (2.12), (2.24) into (2.29), and using (2.25), we obtain T{w,t) in the limit 
r — > oo 



iuis 
8tt ' 



(2.30) 



where we have recalled the definition of s in (1.6). Using (2.8) to extract the Green's 
function from T r we find 

G? 2jl2 (c,0)«-^, (2.31) 

at least in the low frequency limit u — > 0. This is the result claimed in (2.3), giving 
rise to the universal ratio of shear viscosity to entropy density (1.1). 



Acknowledgments 

PB would like to thank the organizers of the ESF Cargese Summer School 2006 "Strings 
and Branes: the Present Paradigm for Gauge Interactions and Cosmology" and the 
organizers of the SNFT 2006 Parma School of Theoretical Physics for the stimulating 
environment created, as well as to all the participants at these schools without whom 
they would have never been the same. AB's research at Perimeter Institute is supported 
in part by the Government of Canada through NSERC and by the Province of Ontario 
through MEDT. AB gratefully acknowledges further support by an NSERC Discovery 
grant. 



9 



A Effective bulk action for ip 

We begin with collecting some useful expressions. The trace of the Einstein equations 
derived from (1.2) takes the form 

R = K a pd<p a d<p b + r ab F$F^ + R+l v . (A.l) 

Additionally we have 



v 7 ^ R x lf ] = - 



c l c 2 C 2 
C 3 



' i , / \(0) 

= ^ V^[V-r ab F^F^j , (A.2) 



where we used superscript (0) to emphasize that the Ricci tensor component and the 
metric determinant are evaluated on the background (1.3). The second equality in 
(A.2) makes use of the Einstein equation for R^ lXl - 

Expanding the effective action (1.2) to quadratic order in perturbation (2.4) and 
using (A.l), we find 

S*akW =T7T^r- [ d D ^y^m \ R m - > 1 ( V - r ab F^F^f } ) 

ib7lLr D+2 Jm d+ 2 {ZD ) 

l — [ d D+ ^V^ {R {2) - ^ R x ^} , 

JM D+ 2 

(A.3) 



16nG 



where in the second line we used identity (A.2). In (A.3), is a perturbation of the 
Ricci scalar, quadratic in tp and its derivatives 

+ ^ R x lf ] , 

where we again used (A.2). From (A.3) and (A. 4) we obtain (2.11). 
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